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The representation of quark distribution and fragmentation functions in terms of non- 
local operators is combined with a simple spectator model. This allows us to estimate these 
functions for the nucleon and the pion ensuring correct crossing and support properties. We 
give estimates for the unpolarized functions as well as for the polarized ones and for subleading 
(higher twist) functions. Furthermore we can study several relations that are consequences of 
Lorentz invariance and of C, P, and T invariance of the strong interactions. 
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I. INTRODUCTION 



Quark distribution functions and quark fragmentation functions appear in the field-theoretical description of 
hard processes as the parts that connect the quark and gluon lines to hadrons in the initial or final state. These 
parts are defined as connected matrix elements of non-local operators built from quark and gluon fields. The 
simplest, but most important ones, are the quark-quark correlation functions [0-|3j. For each quark flavor one 
can write 

$,,(p,p,^) = / e-^-« {p,s\ip^mm\p.s), (1) 

AfcKfc,n,5/0 = E/ (01 {Q\MO\Ph,SH;X) (P,„5,;X|V^,(0)|0), (2) 

diagrammatically represented in Fig. ^ The hadron states are characterized by the momentum and spin vectors 
(limiting ourselves to spinless or spin 1/2 hadrons), |P, 5) and \Ph,Sh)^ for incoming and outgoing hadrons, 
respectively. The quark momenta are denoted p and fc, respectively. In the correlation function A, the sum 
runs over all possible intermediate states that contain the hadron h characterized by Ph and 5*^, while for $ 
the sum is omitted assuming completeness. Furthermore a summation (average) over colors is understood in (f> 
(A). 

In a particular hard process only certain Dirac projections of the correlation functions appear and the non- 
locality is restricted by the integration over quark momenta. For inclusive lepton-hadron scattering, where the 
hard scale Q is set by the spacelike momentum transfer —q^ = Q^, the correlation function $ appears in the 
leading order result in an expansion in 1/Q of the cross section of hard processes. To be precise, the structure 
functions appearing in the cross sections can be expressed in terms of quark distribution functions, e.g., at 
leading order in a^, 2Fi{xb) = F2{xb)/xb = J2a fii^s), where Xb = Q'^/2P ■ q. 

The unpolarized quark distribution function (omitting flavor index a) is obtained from $ as 



I J dp- d^pj. Tt{^j+) 
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FIG. 1. Diagrammatic representation of the correlation functions <1> and A. 
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and depends only on the light-cone momentum fraction x = jf^ / . The lightlike components = a 
defined with the help of lightlike vectors n± satisfying = =^ and n+ ■ n_ = 1. In inclusive lepton-hadron 
scattering they are related to the hadron momentum P and the hard momentum q as 



P = ^ n_ 

Q 



Xb 



V2 



V2 



Q_ 

a/2 



(4) 
(5) 



where M is the mass of the hadron. The non-locality is restricted to a lightlike separation. At this point it 
should be noted that there are other contributions in the leading cross section resulting from soft parts with 
gluon legs. These can be absorbed into the correlation function providing the link operator that renders 
the definition in Eqs. (|l|) and (||) color gauge invariant. Choosing the — gauge in the study of $ the link 
operator reduces to unity. 

The simplest example of a hard process in which the correlation function A appears is 1-particle inclusive 
e^e'" annihilation. In that case the scale is set by the momentum squared of the annihilating leptons, = 
and the production cross section in leading order becomes proportional to fragmentation functions Di^^{zh), 
where Zh — 2Ph ■ q/q^- The quark fragmentation function (omitting flavor index a) is obtained from A as 



Di{z) 



— I dk+ (fkr Tr(A7 
4z j 



Sttz 



k-=p-/z 

Tr{r{^m)\Ph,SH-,X){Ph,Sh-,X\lp{0)\{))) 



(6) 



and depends only on the light-cone momentum fraction z ~ Pf^ /k^. Taking as an explicit example 1-particle 
inclusive electron-positron annihilation, the lightlikc vectors are defined from the hadron momentum P^ and 
the hard momentum q: 



Ph 



ZhQ 

V2 



Ml 



Q 



Zh 



V2 



gV2 



(7) 
(8) 



where is the mass of the produced hadron. The non-locality in the expression for Di is again restricted to 
a lightlike separation. The link operator needed for color gauge invariance becomes invisible by using the gauge 
= in the study of A. 
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As soon as in a hard process two hadrons (or for that matter one jet and one hadron) play a role, the 
transverse directions become important. Examples are inclusive Drell-Yan scattering, 1-particle inclusive lepton- 
hadron scattering or 2-particle inclusive e~^e~ annihilation. For instance, in inclusive lepton-hadron scattering 
the leading order cross section is given by the handbag diagram containing only one soft part <&, shown in 
Fig. I], but in 1-particle inclusive lepton-hadron scattering the leading order cross section involves also the 
fragmentation of the quark, the part A in Fig. 0. In order to describe the current fragmentation one can still 
define lightlike vectors n± using the hadronic momenta as in Eqs. (0) and but the third momentum, in casu 
the hard momentum q, contains a transverse piece. Up to A/^/Q corrections (irrelevant for our purposes), P 
is proportional to is proportional to n_ and 

g=^n_-^rH + ,^. (9) 

By selecting observables depending on the transverse momentum scale —q^ = Qx one needs to consider 

distribution functions and fragmentation functions before integrating over prp or k^, respectively. Examples 
are the dependence on the transverse momentum of the produced lepton pair in Drell-Yan scattering, the 
dependence on the transverse momentum of a produced hadron belonging to the current jet in lepton-hadron 
scattering or the transverse momentum distribution of hadrons with respect to the jet-axis (or with respect to 
a fast hadron in the opposite jet) in the case of back-to-back jets in e^e^ annihilation. 

In this paper, we review the structure of light-cone correlation functions including the effects of transverse 
separation of the quark fields, and we estimate them using a simple model. This will be done for all possible 
Dirac projections that contribute in leading or subleading order. Thus we obtain estimates not only for the 
usual unpolarized distribution and fragmentation functions, but also for the polarized ones and for subleading 
(higher twist) functions. As we will see, in a number of cases there are relations between leading and subleading 
p-jn-integrated functions. 

Although hard cross sections can be expressed in terms of distribution and fragmentation functions, these 
objects cannot be calculated from QCD because they involve the hadronic bound states (at least not for light 
quarks). Even the simplest case, the moments that are related to matrix elements of local operators require 
non-perturbative methods like, for instance, lattice calculations. The full (x-dependent) quark distribution, 
however, requires the knowledge of all moments. We follow here a different route. We want to investigate the 
structure of light-cone correlation functions and illustrate the consequences of various constraints in a simple 
model. Particularly suitable is a model in which the spectrum of intermediate states, which can be inserted 
in Eq. (|l]) or is explicitly present in Eq. (|^), is replaced by one state, referred to as a diquark, if the hadron 
is a baryon. At that point one still has lots of freedom to parametrize the hadron-quark-diquark vertex. We 
make the ansatz that in the zero-binding limit implies the simple symmetric SU (6) spin-flavor structure. As 
parameters to describe the vertex one then has only the quark mass, a diquark mass and a size parameter left. 

Although there are a few parameters, the approach has the advantage of being covariant and producing the 
right support, in contrast to the use of other models B, such as bag models 0-0], quark models |@ M or 



soliton models 11-13|, which require projection techniques ||1J]. Of course, these models have the advantage 
that they also reproduce other observables. Similar spectator models have also been used for the calculation of 
distributions ||l^,^6| and fragmentation functions |l^. 

Of course, matrix elements as discussed above have a scale dependence, which in expressions for hard cross 
sections shows up as a logarithmic dependence on the hard scale. A well known problem is that the modelling 
does not provide a scale dependence. Models with a simple valence quark input, as we take as our starting point, 
will be naturally 'low scale models'. In principle, this could be evolved to higher scales to allow comparison 
with data. However, we lack evolution equations at low scales. On the other hand, we also consider higher twist 
distributions, for which evolution is much more involved than the twist two case JT^ . These are the reasons 
for not including evolution in this paper. In other words, radiative corrections and their absorption in the scale 
dependence of the quark distributions are not taken into account. 

The setup of the paper is the following. In section II we discuss the structure of light-cone correlation 
functions, both quark distribution functions and quark fragmentation functions. In section III we present 
the diquark spectator approach and the construction of the vertices. The results for the distribution and 
fragmentation functions are given in section IV. We end with a summary and outlook. 



II. QUARK CORRELATION FUNCTIONS 
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A. Distribution functions 



In order to study the correlation function $ it is useful to realize that its form is constrained by the hermiticity 
properties of the fields and invariance under parity operation. The most general expression for (f> consistent 
with these constraints is [ p^ , pO[ : 

^p,P,S) = MAi+A2r + A3:^+^at"'P^p, + iA5p-Sj5 + MA6^j5 



(10) 



where the amplitudes Ai depend on cr = 2p ■ P and t = p^ . Hermiticity requires all amplitudes Ai (a, r) to 
be real. Time reversal invariance can also be used and requires the amplitudes A4, A^ and A12 to be purely 
imaginary, hence they vanish. 

In hard processes the hard momentum scale q and the hadron momentum P define the lightlike directions 
n±. The momentum P is parametrized as in Eq. (^. The spin vector S and the quark momentum p are also 
expanded in the lightlike vectors and transverse components: 



Xb{p^+P%) , xQ 

P = H 1= 

xQy/2 Xb\/2 

XMxb XQ 

b = ]=r- n_ ' 



n+ +PT, 



QV2 



MXbV2 



n+ + St- 



(11) 
(12) 



Thus X represents the fraction of the momentum in the + direction carried by the quark inside the hadron. In 
the transverse space the following projectors can be used 

5r=,9'"'-"N-\ (13) 
= e^+A'-. (14) 

Considering a hard scattering process up to order l/Q^, the component of p along n_ is irrelevant and one 
encounters the quantities 

= [[dadrSi )]^'^^^^ 

where we used the shorthand notation 

[dadr 6{ )] = dcrdr S (t - xa + x'^AP + p^) . 

The projections of $ on different Dirac structures define distribution functions. They are related to integrals 
over linear combinations of the amplitudes. The projections 

<^>^-<^\x,PT) = h{x,pl) 

= j [dadrSi )] [A2+XA3] , 



<f[ri {x,pj.)^ i J dp- Tr($r) 



4P+ ' 



(15) 
(16) 



(17) 



[dadrSi )]{ X 



-Ae- 



Pt) 

cT-2x AP 
2M2 



$[^'^'"'^=](x,P7,) = 5^/iiT(a;,p|) 



Pt_ 
M 



Xhi,{x^p?,)+P^-^^ 



M 



(A7 + xAa) 
hiri^^Px) 



M 



(Aj + xAs) 



(18) 



[dadrSi )]{ ^Sir{Ar>+xAw) + 



M 



Ai 



cr - 2a; 
2M2 



An 



p^ prp-Sr 

M M 



An 



(19) 
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are leading in For the distribution functions, this is indicated by the subscript 1 in the names of the 

functions |21[| . The following projections occur with a pre- factor M/P"*", which signals the subleading (or higher 
twist) nature of the corresponding distribution functions 



M 

P+ 
M 

P+ 



[dadrSi )] Ai 



(20) 



$M(^,p^)^^/^(x,p?.) 
M 



AdadrSi )]^^3, 



(21) 



$['''''=l(a;,pj,) 



P+ 
M 

P+ 



P+ 

[dadrSi )] \ -Sip Aq 



M 



M 



S'J~'J)rp iSrpJ)rp 



hrix^Pr) 



— I [dadrSi )] 



M 



2M2 



. Pt Pt-St . 



(22) 



(23) 



$1' 



'(2;,Pt) 



M 



\hL{x,pl) + ^^'f^ hT{x,pir) 



^ j[dadT5{ )]|a 



M 



■A 
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2Af2 

Pt-St fa-2xlVP 



M 



2AP 



a-2x M' 



2\ 2 



2Af2 



An 



An 



(24) 



The constraint of the (5-function in the integration over a and r is indicated in Fig. Furthermore, the 
integration is restricted to the region Af^ = (P — p)^ > 0. This leads to the vanishing of the distribution 
functions at a; = 1. 

If a generic distribution function is written as 



F(x,p|)= [dadrSi )] G {A{a,T), a, x) , 



(25) 



symmetric integration over prp gives 



F{x) = J[dadTe{ )] G{A,{a,T),a, 



where 



[dadr e{ )] = dadr (xa 







(26) 



(27) 



The region covered by the 0- function (for x = 1/2) is the lower shaded region in Fig. g corresponding to 
p'^ > 0. Only the terms involving the distribution functions /i, gi — gi^, hi — hir + {p'^/2M'^) h^rp, e, 
9t = g'x 9t ^^'^ are non-vanishing upon integration over prp. The integrated functions fi{x), 

gi{x) and hi{x) have the well known probabilistic interpretations. fi{x) gives the probability of finding a quark 
with light-cone momentum fraction x in the -I- direction (and any transverse momentum). gi{x) is a chirality 
distribution: in a hadron that is in a positive helicity eigenstate (A = 1), it measures the probability of finding 
a right-handed quark with light-cone momentum fraction x minus the the probability of finding a left-handed 
quark with the same light-cone momentum fraction (and any transverse momentum). hi(x) is a transverse spin 
distribution: in a transversely polarized hadron, it measures the probability of finding quarks with light-cone 
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T = k2 




z = 1/2 / / 
(x = 2)/ / 

A^=i 












o=2k.P 

M^' =(k-P)' 









FIG. 2. The (5-function constraint in tlie cr-r plane (using quark momentum k and hadron momentum P) coming from 
fixing X and fey in the expression for the distribution functions F{x,k'^) (and similarly for the fragmentation functions 
D{z, z'^k'^)) and the full integration regions for the fcr integrated functions F{x) (and similarly for D{z)). The latter 
region is determined by fey > and > 0. 



momentum fraction x polarized along the direction of the polarization of the hadron minus the probability 
of finding quarks with the same light-cone momentum fraction polarized along the direction opposite to the 
polarization of the hadron. The twist three functions have no intuitive partonic interpretation. Nevertheless, 
they are well defined as hadronic matrix elements via Eqs. (|^) and (||), and their projections. 
We note the appearance of higher p|,-moments, 



[dadreO] (^^^-^^—] GiAia,r),a,x), (28) 

such as hf^^'^ and g^^^^ The equality in Eq. ( p8|) is obtained using the azimuthal symmetry of the distribution 
functions, which depend only on x and p^. In the weighted integration, j d^prpptp . . . one will encounter the 

fimctions g[y and h^^\ 

The distribution functions cannot be all independent because their number is larger than the number of 
amplitudes Ai. This is reflected in relations such as 

gT{x)=gi{x) + ^g^^{x), (29) 

hL{x)^h,{x)~^h'[^^\x), (30) 

h^:^\x)^-\^h^^'\xl (31) 

which can be obtained using their explicit expressions in terms of the amplitudes. 
The functions g2 = Qt — gi and /12 = 2(/il — /ii) thus satisfy the sum rules 



L 



dx g2{x)=-g[}^{0), (32) 
dx h2{x) ^2h:l;[^\0), (33) 
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which are a direct consequence of (^) and (po|). If the functions g\rf and h^j^ ' vanish at the origin, we rediscover 
the Burkhardt-Cottingham sum rule and the Burkardt sum rule ||2^. These sum rules (Eqs. ( p2[ ) and ( |33| ) 
with vanishing right-hand sides) can also be derived using Lorentz covariance for the expectation values of local 
operators p3 . In our approach this would imply constraints on the amplitudes Ai . 



B. Fragmentation functions 

The correlation function A is also constrained by the hermiticity properties of the fields and invariance under 
parity operation, leading to an expansion identical to that in Eq. (|lC|) with the replacements {p, P, S', M} 
{k, Ph, Sh, Mh} [0 and with real amplitudes, say Bi, now depending on Th = and ah = 2k- Ph- Time reversal 
invariance does not imply any constraints on the amplitudes, thus B4, and -B12, referred to as 'T-odd', are, 
in general, non- vanishing. (For a discussion on T-odd fragmentation functions, see [p5|). 

In hard processes one encounters the quantities 



A[ri(z,fc^)=l /dfc+Tr(Ar) = f [da^dT^ S{ )] ^ 



^'^^^^ (34) 



where we used the shorthand notation 



[daudTh Si )] = dahdTh S [ - — + ^ + k'j. ] . (35) 



The momentum Ph is parametrized as in Eq. (Q) and is used to define the lightlike vectors, in terms of which 
the vectors Sh and k can also expanded: 



ZhQ z{k^ + k^) 

— 1= n_ H ^— n 1 

zV2 ZhQV2 



(36) 



XhZhQ XhMh , , 

bh = 7= ^ n+ + bhT- (37) 

MhV2 ZhQV2 ^ ^ ' 

Thus z is the fraction of the momentum in the — direction carried by the hadron h originating from the 
fragmentation of the quark. The spin vector satisfies Ph ■ Sh = and for a pure state —S^ = -I- S^rp = 1. 

Comparing the above equations with the case of the distribution functions, one sees that the relations between 
the Dirac projections 2z A^^^{z,kT) and the amplitudes are identical to those for ^^^^{x,prp) after the replace- 
ments {x,a,T,prp,P, St, X,M,Ai,±- components} {1/ z, ah, Th^kr, Ph, ShT, >^h, Mh, Bi,^ - components}, 
except for additional parts originating from the T-odd functions. Furthermore, the definition of fragmentation 
functions follow the general procedure used to define distribution functions. We use for the names of the frag- 
mentation functions capital letters (with the only exception for the coimterparts of /.. functions which are called 
D..). For example, 

A[--](., kr) ^ D,{z, k'S) + '1^2^Di^[z, k'^) 



1 

Yz 



Mh 

[dahdTh (5( )]< 



B2 + -B^ 

z 



^ikr.ShT, (38) 



Mh 



where k'rp — —zkx- The choice of arguments z and kip in the fragmentation functions is worth a comment. In the 
expansion of k in Eq. (|3^ ) the quantities 1/z and kx appear in a natural way. However, in the interpretation 
of A as a decay function of quarks, the variable z as the ratio of Ph is more adequate. Applying a 
Lorentz transformation that leaves the — component (and hence the definition of z) unchanged, one finds that 
k'rp = —zkx is the transverse component of hadron h with respect to the quark momentum. 

Further, we only display the additional parts of projections which come from the time reversal odd amplitudes 
and which, thus, have no counterparts in the distribution functions. In the leading twist projections there is 
only one other projection with an additional term, 

A[-'--=l(z,M--.. + 47^i?i^(-,fe^') 

Mh 

= ... + ^J [dahdTh S{ )] ^ [-B,] . (39) 
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At subleading twist there are five additional T-odd structures: 



P, 



h I . erkTj 



Diiz,k'^) + ei^ShTjDTiz,k'^) 



h ^ 

2zPr 



[datdTh i5( )] < -A,i 



Mh 



Bi2 — e'r^ShTj 



ah - 2Ml/z 



2Ml 



Bi2 



(40) 



A^'''"'"-\z,kT) 



Mh (. p , ,j2, , kT-S, 

— ^XhE,iz,k^) + ^^ 
^ [dahdTh 5{ )] j -A,j 



hT 



2zPr 



ET{z,k'^] 

Gh - 2Ml/z 
2Ml 



krS 



hT 



Mh 



(41) 



^eiH{z,k'^) 
Mh 



2zPr 



[dahdTh 6{ )] < e 



ah - 2Ml/z 
2Ml 



(42) 



The constraint imposed by the 5-function in the ah.-Th plane is also indicated in Fig. The integration is 
restricted to the region = {Ph — fc)^ > 0, which implies that the fragmentation functions vanish dX z = 1. 
We note the reciprocity of x and z, i.e., the constraint for z — 1/2 is the same as one would have for x = 2. 
Note, however, that the integration involves different regions. For the distributions one has (roughly) spacelike 
quark momenta, for the fragmentation timelike quark momenta. If a generic quark fragmentation function is 
given by 



D{z,k'j^) = ^ J[dahdTh6{ )] G{Bi{ah,Th),ah, z), 
the integrated functions are given by 

krp 



(43) 



d'^ki 



2Ml 



TT Z 

"2~ 



[dahdTh e{ )] 



D{z,k'^) 

/ ah-2Ml/z 
V 



G{B^{ah,Th),ah,z), 



where 



[dahdTh 9{ )] = dahdTh T ^ ^'^ 



(44) 



(45) 



Non- vanishing upon integration over fc^ are the fragmentation functions Z^i, Gi ~ Gil, Hi = Hit + 
(fc|/2Af2) H^j,), E, Gt = G'j, + {k^/2Ml) G^, Hl and Dt- 

As for the distributions, the integrated fragmentation functions are not all independent. Using Eq. (^4|) one 
obtains relations such as 



El{z) = z^ 
Dt{z) = 7? 



dz 

d_ 

dz 



e!,'\z) 



D 



-L(l) 
IT 



Griz) = Giiz) - z^ i- 
HLiz) ^ Hi{z) + z^ 



dz 

d 

dz 



g[%) 



H 



z 

-i-(i) 

IL 



- Tz 



H 



(46) 
(47) 
(48) 
(49) 
(50) 
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FIG. 3. The constraint in the a — t plane coming from fixing the spectator mass Mr (compare with Fig. g). 



and similar ones for Dt, G2 = Gt — Gi, H2 = 2{Hl — Hi) and H. Provided that the functions labelled with 
superscript (1) vanish at the origin faster than one power of z, the right hand side vanishes. Finally, let us 
remark that this formalism can be easily extended to include antiquarks pl|]. 



The basic idea of the spectator model is to treat the intermediate states that can be inserted in the definition 
of the correlation function $ in Eq. (|]), or which are explicitly displayed in the definition of the correlation 
function A in Eq. (|^), as a state with a definite mass. In other words, we make a specific ansatz for the 
spectral decomposition of these correlation functions. This may be best illustrated using the support plot 
in (7 and t. In this plot the mass Mu of the remainder, called the spectator, is constant along the lines 
(P — k)^ = T — a + = M^, as indicated in Fig. |3[ The quantum numbers of the intermediate state are those 
determined by the action of the quark field on the state \P, S), hence the name diquark spectator. In the most 
naive picture of the quark structure of the nucleon, such that in its rest frame all quarks are in 1/2+ orbitals, 
the spin of the diquark system can be either (scalar diquark s) or 1 (axial vector diquark a). For a pion state 
we have an antiquark spectator. The inclusion of antiquark and gluon distributions requires a more complex 
spectral decomposition of intermediate states. Here, we restrict ourselves to the simplest case. The correlation 
function $ (the correlation function A will be treated later) is then given in the spectator model by 



project onto different spins in the intermediate state and allow for different spectator masses. 

We start with the correlation function $ for a nucleon. The matrix element appearing in the RHS of ( |5^ ) is 
given by 



leading to 




(51) 



III. THE SPECTATOR MODEL 





(53) 
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in the case of a scalar diquark, or by 

(X(^)|^.(0)|P,5) =e;(^)(^^y^ T^f UdP,S), (54) 

for a vector diquark. The matrix elements consist of a nucleon-quark-diquark vertex T(iV) yet to be specified, 
the Dirac spinor for the nucleon Ui{P, S), a quark propagator for the untruncated quark line (m is the constituent 

mass of the quark) and a polarization vector e^^^^ in the case of an axial vector diquark. The next step is to fix 
the Dirac structure of the nucleon-quark-diquark vertex T. We assume the following structures: 

T^(iV) - 1 5s(p'), (55) 



The functions gnip^) (where i? is s or a) are form factors that take into account the composite structure of the 
nucleon and the diquark spectator. In the choice of vertices, the factors and projection operators are chosen to 
assure that in the target rest frame, where the nucleon spinors have only upper components, the diquark spin 
1 states are purely spatial and in which case the axial vector diquark vertex reduces to xlv"" ' ^Xq- The most 
general structure of the vertices can be found in [E6[ . With our choices, we find 



P, S) = \^ ^^^'7^;/^^^ (1^ + ™) (f + M) (1 ^ a.,,^) + „.), (57) 



where is a spin factor, which takes the values = 1 and aa — —1/3. In obtaining this result we used as the 

polarization sum for the axial vector diquark in the form e^'"'^''e!>^' ~ ^gfif + P^i.Pv/M'^ , which is consistent 
with the choice that the axial vector diquark spin states are purely spatial in the nucleon rest frame. We will 
use the same form factors for scalar and axial vector diquark: 

9ir)-N^-^. (58) 

The quantity A is another parameter of the model which ensures that the vertex is cut off if the virtuality of the 
quark leg is much larger than A^. A'^ is a normalization constant. This choice of form factor has the advantage 
of killing the pole of the quark propagator as suggested in . 

In the same way, one can write down a simple spectator model for the pion. The matrix element can be 
written as 

(X(")|V,(0)|P,> = f-^) Tkiv[''\ (59) 

The spinor v'f^ describes the spin state of the antiquark spectator. The simplest vertex is given by 

T(.) = ^^±^75. (60) 

Taking for the spectator antiquark spin sum ^^1'"''^/"'' = /jt ~ M-,^ one arrives at precisely the same expression 
as for the nucleon (Eq. (p7|)) with — Q. 

From the correlation function $ one easily obtains the distribution functions. Taking out the explicit 5- 
function, 

<i>«(p, P, S) = $(p, P,S)5{T-a + A/2 - Ml) , (61) 
one finds immediately from Eq. (|l^) the result 



^rr,, , Tr($r) 



(62) 



with 
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Pt 



1 



1 



■Mi 



(63) 



We now turn to the fragmentation functions. The calculation is very similar to the case of the distribution 
functions, involving the same type of matrix elements. Further, we assume that the hadron h has no interactions 
with the the spectator. This allows us to use a free spinor to describe this outgoing hadron. Then we see that 
the correlation function A is the same as the one needed for the distributions, after obvious replacements in the 
arguments, namely 



2(27r)3 



+ m) ifh + Mh) (1 + a^jrjf,) i^ + m). 



(64) 



A direct consequence is 



Ari(z,M = ;^<I'["'- 
2z 



1. $[r'] 
2z 



(65) 



where V and F involve an interchange of 
A(fc, Ph,Sh) S {{k - Phf - Ml), Eq. dH) leads to 



and 



components. Writing A{k, Ph, Sh) 



A[ri(z,feT) 



Tr(AF) 



(66) 



with 



1- z 



Mi 



(67) 



The consequence of using free spinors to describe the outgoing hadron is that all T-odd fragmentation functions 
vanish and we have a one-to-one correspondence between distribution and fragmentation functions. As can be 
seen in Fig. ^ the actual behavior of the distribution and fragmentation functions comes from different regions 
in T, roughly spacelike and timelike, respectively. Therefore, the above reciprocity (Eq. (|65[)) is of use for the 
analytic expressions, less for the actual values. 



11 



IV. RESULTS AND DISCUSSION 



A. Distribution functions of the nucleon 



Using the expression in Eq. { p7\) we can compute the amphtudes Ai shown in Eq. (10) 
common factors by defining 



S{t-(j + M2 - Ml) 



:2a 



A,. 



2(27r)3 |r-A2| 
we obtain, as expected, the T-odd amphtudes A4 = A5 = A12 = 0, and 



Ai 


^ ((M + mf - 


Ml) 




A2 


= -{r- 










= {M + 


mf - Ml 


+ {r- 








VM 




Ml) - 


A7 


= 2 Ah mM, 






is 








ig 


= an (r 








Aw 




\M + mf 


-Ml 




All 


= -2aRM^. 







m 
M 



m 
M 



Introducing the function A^(a;) such that 



,2 2_PT + ^li^) 



1 -X 



which imphes 

Xl(x) = A2(1 ~x)+ xMI - x{l - x)M'^, 
one gets the following results for the distribution functions, 
iV^ (1 _ a;)2a-i [xM + mf+p^j, 



fi{x,PT) 



9il{x,Pt) = an 



16-^ {Pl + Mf- ' 

iV2(l-x)2«-i {xM + mf~p^^ 
16^r^ [pl + \l) 

(1 - M{xM + m) 
8^^^ (p| + A2/" 



2a 



hlTix,PT) = aa fi{x,p^), 



hhix^Pr) = -girix^Pr), 



Kt{x,Pt) ^ -0.R 
e{x,p'^) 



(1 - x) 



2a-l 



M2 



8.3 (p| + A|/"' 

7V2 (1 _ x)2"-2 (1 - a;)(a;Af + m)(Af + m) - Ml (a; + f ) - (l + f ) p^. 



Taking out some 



(68) 



s2a 



^2 _ .x)2"-2 (1 _ a;2)M2 + 2mM(l - x) - Af|, - p^ 



16.3 



(69) 

(70) 
(71) 
(72) 

(73) 
(74) 
(75) 
(76) 
(77) 



(78) 
(79) 

(80) 
(81) 

(82) 
(83) 

(84) 
(85) 
(86) 
(87) 
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(88) 



/i^(x,p|) = a«/^(x,p|,), (91) 
2 jV^ (1 - :r)^°-^ (1 - x){xM + m)(M + m) - {x + §) + (l - 2x - f ) 

hT{x,p?r) = -9i{x,pl). (93) 

Although there is a certain freedom in the choice of the parameters, one immediately sees that the occurence 
of singularities in the integration region (see Fig. ^ will cause problems which are avoided if there is no zero in 
the denominator. The requirement that A^(a;) is positive implies for the distribution functions (0 < a; < 1) 

Mr > M-A, (94) 
while for the fragmentation functions (using reciprocity, we have to look at a: > 1) it leads to 

Mii>A-Mh. (95) 
Provided condition (M) is fulfilled, one obtains the integrated distribution functions, 



jy2 _ a;)2a-i 2(a - 1) (xM + m)2 + Xl{x) 
32^2 (a_i)(2a-l) (A2(a;))' 



_ N^a^(\~xf^-^ 2{a~\){xM-rmf~\l{x) 
^'^''> 32^2(«-l)(2«-l) (AKx))^"-^ ' ^ ' 



N^a, (1 - x)2"-i {xM + m)2 
"^^^"^ - 16.^(2a-l) (Ai(.))--- ^''^ 



_ (1 ~ a;)2"-2 2(a -l){x+ t) [(1 - x)(M + m)Af - M^] - (l + f) A|(x) 

"^"^^ - 32^2(^ _ 1)(2„ _ 1) (A|(x))^"-^ ' ^ ' 

jW (l-a;)2°-2 2(a - 1) (x + t) [(1 ~ x)(Af + m)M - M^] - (x + f ) Ai(x) 

^^^^^ ^ 32.>-l)(2«-l) ^^T^^jf^-^ ' (100) 

, , , _ A^^flH (1 - xr^-^ 2{a -l){x+ t) [(1 -x){M + m)M- M|] + (l - 2a: - f ) Ai(x) 

'^^^^^ - 32.2(, „ i)(2„ _ 1) {Xl{x)f''-' ■ ^ ' 

Examples of the p|./2M2-weighted distributions are 



327r2(a-l)(2a-l) ^XUx))' 



We note that these latter functions do not vanish at a; = 0, implying non- vanishing sum rules for 52 and in 
accordance with Eqs. ( ^2|) and (^), except if the quarks are massless. 
The functions g2 and /12 are given by 

h,{x) N^a, (1 - a:)2°-2 2{a - I) {x + ^) [M\1 - xf - M|] + [l ^ 2x - fj) Xjjx) 

2 32^2(„_i)(2«_i) (A2(:.))^«-i ' ^ ' 



We can directly check that Eqs. ( |29| ) and (30) are satisfied. 

Up to now, we have not specified flavor in the distributions. For the nucleon we only distinguished two types 
of distributions, fl and /°, etc. Since spin diquarks are in a flavor singlet state and spin 1 diquarks are in 
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a flavor triplet state, in order to combine to a symmetric spin-flavor wave function as demanded by the Pauli 
principle, the proton wave function has the well-known SU{4:) structure. 



b T> = ^ I" T 5°) + -^\u t A") - i |y i Al) - i Id t A?) + Id i Al), 



(104) 

where S (A) represents a scalar (axial vector) diquark and the upper (lower) indices represent the projections 
of the spin (isospin) along a definite direction. Since the coupling of the spin has already been included in the 
vertices, we need the flavor coupling 



\P) = ^ I" + ^o) - ^ M ^i>, 

to find that for the nucleon the fiavor distributions are 

S 1 
Ji ^ 2 1 2 

rd j'a 

Jl — Jl : 



(105) 

(106) 
(107) 



and similarly for the other functions. The proportionality of the numbers is obtained from Eq. ( |105| ), while 
the overall factor is chosen to reproduce the sum rules for the number of up and down quarks if fl and are 
normalized to unity upon integration over p-^ and x. This will fix the normalization N in the form factor in 
Eq. (|5|). Notice that the factors — I and aa = —1/3 in the distribution functions will produce different u and 
d weighting for unpolarized and polarized distributions. Further differences between u and d distributions can 
also be induced by different choices of M^, A or a. We take for the nucleon a = 2 to reproduce the right large x 
behavior of i.e. (1 — x)^, as predicted by the Drell-Yan-West relation and reasonably well confirmed by data. 
We refrain from tuning the large x behavior of ff to match the (1 — x)'^ form indicated by data. Since is 
only affected by vector diquarks, this could be easily obtained by choosing a different form factor for the latter. 
We feel that this kind of fine-tuning would take things too far with the simple model we use. Similarly, we will 
only consider one common value for A. We will, however, consider different masses for scalar and vector diquark 
spectators. The color magnetic hyperfine interaction, held responsible for the nucleon-delta mass difference of 
300 MeV, will also produce a mass difference between singlet and triplet diquark states. Neglecting dynamical 
effects, group-theoretical factors lead to a difference Ma — Ms — 200 MeV p^ . 

Another important constraint comes from the axial charge of the nucleon, given by 



9a = [\x[gnx)-9fix)] = 

^0 "'0 



dx 



(108) 



The sensitivity to the parameters Mji and A is best illustrated by considering some characteristic values. We 
take a quark mass of 0.36 GeV (about one third of the average nucleon-delta mass), two different values for M/j 
(0.6 and 0.8 GeV) and three values for A (0.4, 0.5 and 0.6 GeV). The distributions turn out to be insensitive to 
the value chosen for the quark mass. In Table | the values of some moments are given. 

Fig. ^ shows the twist two distributions for different values of the mass of the spectator and of the parameter 
A. Clearly, AIr dictates the position of the maximum, while A governs the width of the distribution. We can 
see that an increase of Mr induces a shift on the peak of the valence distribution fi{x) towards lower values of 
X and a decrease in its second moment. In order to model sea quark distributions, one could introduce heavier 
spectators with four quarks or three quarks and one antiquark, and in this way satisfy the momentum sum rule at 
the model-scale. In the absence of transverse momentum for the quarks he have fi{x) = gi(x)/aji — hi{x)/ a^. 



TABLE I. The second moment of /i, {x) = J dxxfi{x) and the first moments gi ~ J dx gi{x) and hi — J dxhi{x) 
are given for two diquark masses and for three values of A. 



A (GeV) 


(^>^ 


Mr = 0.6 GeV 

9? 




{x)^ 


Mr = 0.8 GeV 




0.4 


0.366 


0.923 flfi 


0.962 Ofl 


0.230 


0.650 aR 


0.825 OH 


0.5 


0.375 


0.794 flH 


0.897 Ofl 


0.256 


0.527 ajj 


0.764 Ofl 


0.6 


0.384 


0.671 an 


0.835 Ofl 


0.277 


0.416 ajj 


0.708 Ofl 
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TABLE II. The second moment of /i, {x) = J dx xfi{x) and the first moments gi — J dx gi{x) and hi = J dxhi{x) 
are given for it and d quarks in a proton for three values of A. 



A (GeV) 




gi 




{xf 




hi 


0.4 


0.664 


1.277 


1.305 


0.230 


-0.217 


-0.275 


0.5 


0.690 


1.103 


1.218 


0.256 


-0.176 


-0.255 


0.6 


0.715 


0.937 


1.135 


0.277 


-0.139 


-0.236 



In Table y we have given the values of some moments for u and d quarks in the proton using Mg = 0.6 GeV 
and Ma = 0.8 GeV. We now use the axial charge of the nucleon, gA — 5" — gf, to find the most suitable values 
for A. The value A = 0.5 GeV gives g^ = 1.28, close to the experimental value. 

Fig. H shows the distribution fi{x) multiplied by x. We find a satisfactory qualitative agreement with the 
valence distributions of Gliick, Reya and Vogt (GRV) calculated at the low scale ^j^Q = 0.23 GeV^ ||2^. For 
u and d quarks, the first moment of fi{x) is clearly larger in our model, which would imply that our results 
describe the nucleon at an even lower scale than GRV. Fig. ^ shows the distributions gi{x) and gf{x) multiplied 
by X for the values Ms = 0.6 GeV, Ma = 0.8 GeV and A — 0.5 GeV. Again, we find a qualitative agreement 
with the polarized valence distributions of Gliick et al. |2^. Using the same parameters, we can obtain higher 
twist distributions. Fig. |^ shows the twist three distributions e"(a:) and e'*(x), Fig. ^ shows the distributions 
g'^^^^{x) and g''y^'^{x)^ while in Fig. ^ we have the distributions 52 (^) qH^)- ^^'^^ * small violation of 
the Burkhardt-Cottingham sum rule, in agreement with Eq. (|3^), due to the non-zero value of g^^{Q). 

At this point it is important to realize that in this model there are no antiquarks. This means that the 
distribution functions are zero for x < 0, due to the symmetry properties of the matrix elements involved in the 
calculation. Therefore, C-even and C-odd sum rules are equal. 



B. Distribution functions for the pion 



The expressions for the distribution functions for the pion are the same as those for the nucleon with 0^ = 
and making the replacement Mji — > m. Only spin independent functions will remain. For a = 1 we have fi{x) = 
2(1 — x), showing the power law behavior expected from simple counting rules. The apparent singularity caused 
by the factor (a — 1) that enters in the denominator can be avoided including this factor in the normalization 
N . In this case we find the symmetry x ^ (\ ~ x) for xfi jx). For values of a different from 1 the functions 
depend on the parameter A, which is constrained by Eq. (Q). For A — 0.4 GeV, the distribution fi{x) is shown 
in Fig. ^ for two values of a and compared with the parametrization of the leading order valence distributions 
of Gliick, Reya and Vogt at the low scale — 0.25 GeV^ 

In this case the vertex gives immediately identical antiquark distribution or, equivalently, a contribution for 
negatives values of x. 



C. Fragmentation functions for the nucleon 



The assumed form of the quark-diquark-nucleon vertex also allows the calculation of the fragmentation 
functions for the nucleons. We can use the reciprocity relation mentioned at the end of section III to obtain the 
analytic expressions for Di(z, —zkx), etc., and, after integration over krp = —zkx, the expressions for Di{z) 
etc. In this case we introduce the function Ah through the relation 

For example, the unpolarized fragmentation function reads 

L)i[Z,ZkT)- ^g^3^2a — .„ .ixx2a' i^-^Uj 
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TABLE III. The first two moments of the unpolarized fragmentation function {N^ = J D^{z)dz and e^) and 
the ratios G^/N^ = / Gf (2) dz/N^ and H^/N^ = J H^{z) dz/N" for different values of A and Mr. 







Mr = 


0.6 GeV 








0.8 GeV 




A (GeV) 




e« 




Gf/iV« 




Nr 










0.4 


1.908 e" 


2.732 X 10~ 


5 


0.628 ur 


0.814 flR 


2.112 e"" 


3.785 X 10" 


4 


0.538 aR 


0.769 Oh 


0.5 


1.901 


4.808 X 10" 


4 


0.639 ur 


0.820 flR 


2.107 e'^ 


1.345 X 10" 


3 


0.548 aR 


0.774 Oh 


0.6 


1.891 e"" 


2.429 X 10" 


3 


0.654 aR 


0.827 aR 


2.099 e"^ 


3.732 X 10" 


3 


0.561 an 


0.781 Oh 



and, after integration over the transverse momentum, 

Nl (1 - 2{a - 1) {M + mzf + (1) 



327r2 (a - l)(2a- 1) 



(-^Ai(i)) 



2q-1 



(111) 



The factor Np is a normahzation constant. Distinguishing D\(D'{) as the fragmentation functions for a quark 
into a nucleon and an anti-S diquark (anti-A diquark), one finds 



^1 J 



(112) 
(113) 



and similarly for G\ and Hi. In this case there is no sum rule to fix the normahzations of D\ and D\. In the 
5t/(4) symmetric limit (D\ = D^), one finds the expected ratio D"^^ / D^^^ = 2. We will introduce the scale 
invariant quantities 



dzzDi'iz), 



(114) 



and express our results with the help of these quantities. In Table EII the values of some moments are given for 
a quark mass of 0.36 GcV, two different values of (0.6 and 0.8 GeV) and three values of A (0.4, 0.5 and 0.6 
GeV.) By normalizing e** = e° = 1, we obtain the results for u and d quarks given in Table IV. 

In Fig. |ll| we show our results for the unpolarized fragmentation function Di and the ratios of polarized to 
unpolarized functions Gi{z)/ Di(z) and Hi(z)/ Di(z). Since the dependence on the parameter A turns out to be 
weak, we display results for the choice A = 0.5 GeV only. On the left hand side of the figure the two spectator 
masses {Mr = 0.6 GeV and Mr = 0.8 GeV) are compared. On the right hand side we show the results for 
the u and the d quark fragmentation functions as defined by Eqs. (112) and ( |113| ). To allow a comparison with 
data |31| we fixed the normalization such that the second moment e?"*?" — J dz z D\{z) takes the value 



dzz (D'i^Piz) - D'i^f'iz)) « 0.019, 



(115) 



which is our (rough) estimate for the second moment obtained from the EMC data. We compare our result for 
Di{z) to the difference 



Z?P^'(z) - i?pP(z) = i?pP(0) 



(116) 



TABLE IV. The first two moments of the u and d quark combinations of Eqs. (112) and (113) and the ratios 
Gl/N" = / Gl{z)dz/m and Hl/N" = J Hf{z)dz/N'' for fragmentation into protons. 



A (GeV) 


iV" 




u-quark 

Gr/A'" 


Hi" /A" 






d-quark 

Gf/N'' 


Hf/N'^ 


0.4 


3.917 


2 


0.410 


0.525 


2.112 


1 


-0.179 


-0.256 


0.5 


3.904 


2 


0.418 


0.529 


2.107 


1 


-0.183 


-0.258 


0.6 


3.886 


2 


0.427 


0.533 


2.099 


1 


-0.187 


-0.260 
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since it is the appropriate combination for comparison with a model involving only valence quarks, although in 
our model Z)"~*^ is zero by construction (as are all so-called unfavored fragmentation functions). Furthermore, 
since in the EMC analysis the assumption D^^^{z) — D'l^^{z) was made, we compare our result for Df{z) 
to just the same combination of Eq. (116). The ratios of polarized to unpolarized fragmentation functions, 
Gi{z) / Di{z) and Hi{z) / Di{z), are given for the q ^ u,d as well; normalization factors drop from the ratios. 



D. Fragmentation functions for the pion 

Our results for the fragmentation function of a quark to a pion are shown in Fig. |l^. We display zD^^'^* (z) 
which is rescaled such that the second moment e'^^'" = f dz z D^^^* (z) equals the value of the valence combi- 
nation 

J dzz (Di'^'^^ (z) - D'^^''' (z)) « 0.088, (117) 

our estimate for the second moment obtained from the corresponding EMC data pl| | (more recent parametriza- 
tions available for the combination {D'^^''^ + D'^^''') ^ agree with the EMC data). 
Note that in our calculation all favored fragmentation functions are identical 

D'^^^'^iz) = Df-^''*{z) = Df-^'^'iz) = Df^'^'iz), (118) 

while all unfavored fragmentation functions have not been considered in our approach: 

Df-^""^ {z) = Df^""* (z) = Df-^""' (z) = D'^^''' (z) = 0. (119) 

The latter property has to be contrasted with the experimental observation that unfavored fragmentation 
functions can be as large as the favored ones for small z, are suppressed by a factor of about 2 for 0.4 < z < 0.6 
and even stronger suppressed for large z [pi[. This observation holds true for both nucleons and pions. 



E. Summary 

In this paper we combined the representation of distribution and fragmentation functions in terms of non-local 
operator expectation values with a simple spectator model. This amounts to saturating the antiquark-hadron 
intermediate state with one single state of definite mass. With an effective vertex that connects to this state, 
containing a form factor, we can calculate all the independent amplitudes that are allowed for the non-local 
operator expectation values after imposing constraints of Lorentz invariance and invariance under parity and 
time-reversal operations. Exploiting the explicit expressions of distribution and fragmentation functions in 
terms of those amplitudes, several relations between p^-integrated distribution (or fragmentation) functions 
arise. 

For nucleons and pions we have obtained expressions for the distribution and fragmentation functions within 
our approach. Flavor charges and axial vector charge served to fix the free parameters of the model in the 
case of the distribution functions. For the fragmentation functions we utilized the same set of parameters 
except for the overall normalization which in this case is not constrained by a number sum-rule. Considering 
all Dirac projections in leading and subleading order (in an expansion in 1/Q) we were able to give estimates 
for the polarized and unpolarized cases, including subleading (higher twist) functions. The latter lack a simple 
probabilistic interpretation, but are well defined as projections of non-local operators. 

By comparing our expressions to available parametrizations at 'low (hadronic) scales' and to some experimen- 
tal data we find that we can obtain reasonable qualitative agreement for the unpolarized distribution function 
/i(x), the longitudinal spin-distribution 51(2^) and with the unpolarized fragmentation function -Di(z) for both 
nucleons and pions. These findings give confidence that the estimates obtained for the 'terra incognita' func- 
tions (transverse spin distributions, longitudinal and transverse spin fragmentations and subleading functions) 
provide a reasonable estimate of the order of magnitude of the functions and their (large) x behavior despite 
the simplicity of the model. The comparison to the available parametrizations and experimental data gives an 
indication of the level of accuracy our estimates can reach, keeping in mind that we have excluded the sea-quark 
and gluon sectors, and evolution. 
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FIG. 4. Twist two distributions for the nuclcon. The plots on the top represent fi{x), the ones on the middle show 
gi(x)/aR and at the bottom we have hi{x)/aR. The plots on the left correspond to A = 0.4 GeV and the ones on the 
right to A = 0.6 GeV. The full line corresponds to Mr = 0.6 GeV and the dashed line to Mr = 0.8 GeV. 
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0.5 1 




0.5 1 




0.5 1 



FIG. 5. Twist two distributions for the nucleon. The plot at the top shows xfl{x) (full line) and xfi{x) (dashed line) 
for Ms = 0.6 GeV, Ma = 0.8 GeV and A = 0.5 GeV. The plot on the middle shows xfi'ix) (full line) and xft{x) (dashed 
line) for the same values of the parameters. The third plot shows the low scale {p'^ = 0.23 GeV^) valence distributions 
of Ghick, Reya and Vogt 
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Xgi(x) A 




-0.5 



Xgi(x) A 




-0.2 - 



FIG. 6. Polarized proton distributions gi{x) and gtix). The first plot shows our estimates for xgx{x) (full line) and 
xgi{x) (dashed line) for A = 0.5 GeV, Ms = 0.6 GeV and Ma = 0.8 GeV. The second plot shows the low scale = 0.23 
GeV^ parametrization of Gliick, et al. Eol for the same functions 
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FIG. 7. Distributions e"(a;) (solid line) and e''{x) (dotted line) for Ms = 0.6 GeV, Ma = 0.8 GeV and A = 0.5 GeV. 
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82M 




FIG. 9. Distributions gS{x) (solid line) and gi{x) (dashed line) for = 0.6 GeV, Ma = 0.8 GeV and A = 0.5 GeV. 
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FI G. 1 1. Twist two fragmentation functions. The plots on the left show D^{z)/e^' (with defined according to 



Eq. ( |114D ) and the ratios Gi(z)/{aRDi^{z)) and Hi(z)/(aflDf (2)) for Mr = 0.6 GeV and Mr = 0.8 GeV (A = 0.5 
GeV). The plots on the right show the u and d quark results obtained with = e° = 1. For comparison with data our 
results for Dl(z) are rescaled (see text). Data for the difference (D"~*^(2) — D^~'^{z)) are taken from |3l| ]. 
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